Possible routes for synthesis of new boron-rich Fe-B and Fei_ x .Cr x .B 4 compounds 
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We use ab initio calculations to examine thermodynamic factors that could promote the formation 
of recently proposed unique oP10-FeB4 and oP12-FeB2 compounds. We demonstrate that these 
compact boron-rich phases are stabilized further under pressure. We also show that chromium 
tetraboride is more stable in the new oPIO rather than the reported oIlO structure which opens up 
the possibility of realizing an oP10-(Fe I Cri_ :r )B4 pseudobinary material. In addition to exhibiting 
remarkable electronic features, oP10-FeB4 and oP12-FeB2 are expected to be harder than the known 
Fe-B compounds commonly used for hard coating applications. 



The renewed interest in transition metal (TM) borides 
stems, in part, from the materials' potential to serve 
as hard, wear-resistant, chemically inert coatings^— 
Combination of the TM and boron ensures a high va- 
lence electron density and a pronounced covalent bond- 
ing resulting in the compounds' exceptional hardness; 
for example, ReB2 has been recently demonstrated to 
be the first metal-based bulk material that can scratch 
diamond*^ Metallicity and covalency are also key ingre- 
dients for phonon-mediated superconductivity and TM 
borides have received a lot of attention^ following the 
discovery of a remarkable MgB2 superconductor iL 

Iron borides have two particularly important industrial 
applications. First, homogeneously dispersed second- 
phase particles of Fe2B are known to improve the ten- 
sile strength of low-carbon steels^ while pseudobinary 
Fe3(B,C) and Fe23(B,C)g2»i£ precipitates harden high- 
carbon steelsJ^ Second, Fe-B-based hard protective coat- 
ings are produced directly on the surface of steel via 
the process of boridingJ^— During this thermochemi- 
cal process boron diffuses into the steel forming either a 
single-phase (Fe 2 B) or a duplex-phase (Fe2B+FeB) coat- 
ing layer The two FeB and Fe2B compounds have been 
shown to crystallize in the oP8 (or the related oS8)^ 
and tI12 configurations, respectively, and are the only 
low-temperature ground states listed in the latest ex- 
perimental phase diagram^ Synthesis of new boron-rich 
Fe-B compounds could have technological implications 
as the hardness of metal borides tends to increase with 
boron content. So far, except for the observation of 
a metastable FeB 4 g intercalation compound^ only two 
studies reported on the synthesis of amorphous^ and 
the AlB2-typei^ iron diboride but these compounds have 
not been reproduced. 

Application of advanced compound prediction meth- 
ods has recently allowed us to identify oP12-FcB2 and 
0PIO-FCB4 candidate ground states with unique crys- 
tal structures (Fig. lb,c) that are stable relative to the 
known compounds^ oP12-FeB2 was predicted to be the 
first metal diboride semiconductor while oP10-FeB4 was 
shown to have the necessary features to exhibit phonon- 
mediated superconductivity with a T c of 15-20 K. The 
aim of this study is to examine synthesis routes that 
could lead to the discovery of the new materials. We 



find that both compounds arc stabilized further under 
pressure while the oPIO Fe-based phase could be real- 
ized in the (Fe a; Cri_ :E )B4 pseudobinary form. We also 
observe that calculated elastic constants of the predicted 
Fe-B compounds are higher than those in the known FeB 
and Fe2B materials. 

We use the projector augmented wave method 2 ^ as im- 
plemented in VASR2 2 - and carry out full structural and 
spin relaxation for all compounds; relevant A/B2 and 
M B4 metal borides show no magnetic ordering and their 
phonon and electron-phonon (e-ph) properties arc exam- 
ined without spin polarization. The chosen energy cutoff 
of 500 eV and dense Monkhorst-Pack fc-meshes 2 ^ ensure 
numerical convergence of formation energy differences to 
typically 1-2 meV/atom. We employ the Perdew-Burkc- 
Ernzerhof (PBE) exchange-correlation (xc) functional 2 ^ 
within the generalized gradient approximation (GGA). 
Vibrational corrections to Gibbs energy are calculated 
with PHONj2£ The strength of the e-ph coupling for 
oP10-CrB4 is evaluated within the linear response the- 
ory using the Quantum-ESPRESSO packaged The boron 
ground state is simulated as a-B, known to be stable 
at medium pressures up to ~ 19 GPai 27 ' 28 The iron 
ground state is taken as bec and hep at and 20 GPa, 
respectively^ 

The peculiar 0IIO-TMB4 phases comprised of tetrag- 
onal B nets (Fig. la) were previously examined within 
the extended Hiickel method; it was concluded that max- 
imum binding in the 3d scries is achieved for Cr and that 
the electron-rich Fc, Co, and Ni tctraborides may be un- 
stable in this configuration^ Our calculations showed 2 ^ 
dynamical instability of the oI10-FeB4 phase which gains 
a considerable 0.13 eV/f.u. in enthalpy by transforming 
into the oPIO structure (Fig. 2b). In the present study 
we observe a similar behaviour in the Cr-B system: oIlO- 
CrB 4 is found to be both dynamically (Fig. 1 in Ref. l3ll ) 
and thermodynamically (by 0.03 eV/f.u.) unstable rela- 
tive to oPIO. Although the boron network undergoes a 
significant distortion in the oIlO— >oP10 transformation, 
the lattice parameters and the simulated powder diffrac- 
tion patterns (Fig. 2 in Ref. l3lh remain close which 
may explain why the structure of &B4 was originally 
solved as oIlO. These findings suggest that Fe^Cri-j^ 
compositions may assume the oPIO structure as well un- 
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FIG. 1: (Color online) a-c) Boron-rich structures: iron and 
boron are shown as large cyan and small black spheres, respec- 
tively, d) Atomic volume as a function of pressure, e) Relative 
enthalpy of oP12-FeB2 and oP10-FeB4 candidate phases w.r.t. 
the a-B<r->-oP8-FeB tie-line. /) Relative Gibbs energy (with 
vibrational contributions included using PHON^) of oPlO- 
FeB4 w.r.t. a-B<->oP12-FeB 2 at K with and without zero 
point energy (ZPE) and at 1000 K. 



der standard synthesis conditions. Indeed, there were 
no special requirements regarding the heat treatment or 
the starting materials (apart from cold compacting of the 
MB powder s 32 ' 33 ) for the past synthesis of ternary Fe- 
Cr-B materials, e.g., the ordered metal-rich Mn4B-type 
alloy^i or the oP8-FeB- and oS8-CrB-type compounds 
with partial substitutions of the host metah 32 ' 33 

The effect of the composition on the thermody- 
namic and electronic properties of the ordered oPlO- 
Fe :r Cri_; E B4 compounds is investigated via supercell 
simulations^ see Fig. 2. We observe a sizeable sta- 
bilization, up to 28 meV/f.u., of the oPIO structure at 
x = 0.5, Fig. 2a. The configurational entropy contri- 
bution from the disorder in population of metal sites 
is comparable at elevated temperatures: AG cora /(T) = 
— kBT[x\nx + (1 — x) ln(l — x)} would be 60 meV/f.u. at 
T=1000 K and x = 0.5 if all decorations were degenerate 
in energy. Hence, the resulting ordering of Fe and Cr on 
the metal sublattice will depend on the quenching condi- 
tions. The density of states (DOS) is found to be sensitive 
to metal site population patterns at the Fe-rich end and 
drops rapidly in going from pure FeB4 (1.0 states/(spin 
cV atom)) to pure CrB4 (0.18 states/(spin eV atom)), 
Fig. 2b. Our linear response theory calculation s 26 ' 36 give 
a small e-ph coupling (A « 0.15) and a negligible critical 
temperature (T c C 1 K) in oP10-CrB4 indicating that 
superconductivity in this pseudobinary will require high 
Fe concentrations. 

Application of medium pressures (a few GPa) in multi- 
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FIG. 2: (Color online) a) Relative stability of the pseudobi- 
nary oP10-Fe I Cri_ I B4 compound, b) Total and projected 
density of states (DOS) for lowest enthalpy compounds. 



anvil or diamond anvil cell setups may promote the for- 
mation of the materials in the predicted configurations by 
improving the compounds' thermodynamic stability or 
the reaction kinetics. Successful examples of this synthe- 
sis route include boron-rich CaB4 and NdBg or metal-rich 
Fe2B materialsi 3 -^— We have examined the response to 
the hydrostatic pressure of over 40 known and proposed 
ambient-pressure M-B structure types listed in Ref. Hi] 
by calculating their formation enthalpies at P = 20 GPa. 
We find that oP8-FcB is the lowest-enthalpy phase and 
that the necessary condition for a boron-rich phase to 
be thermodynamically stable is to lie below the 0P8- 
FeB-H-cv-B tie-line. Due to the remarkable compactness 
(Fig. Id) of the predicted oP10-FeB 4 and oP12-FeB 2 
phases (a respective 5% and 7% reduction in atomic vol- 
ume compared to mixture of a-B and FeB) pressure does 
lead to much lower relative formation enthalpies for both 
compounds (Fig. le). Above P = 10 GPa, oP10-FeB 4 
becomes thermodynamically stable relative to a-B and 
oP12-FcB2 at all temperatures (Fig. If). 

The compactness of these Fe-B phases also opens the 
perspective of obtaining new metal-based hard materi- 
als. It has been argued that the extraordinary hardness 
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TABLE I: Calculated bulk (B), shear (G), and Young (E) 
moduli and Poisson ratio (v) for known and predicted Fe-B 
phases. hP6-ReB2 is added for comparison. 
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of the ReB2 material arises from an efficient packing of 
B into hep-Re that results in only a 5% expansion of 
the metal lattice and, consequently, in the shortest TM- 
TM distances for any TM diboride»i The importance 
of having strong TM-TM and TM-B bonds is that the 
(0001) plane in hP6-ReB2 supports the weakest stress. 
In the predicted oP12-FeB 2 , the TM-TM (TM-B) bonds 
are shorter by 11% (9%) not only because of the smaller 
size of Fe but also because of the break-up of the B lay- 
ers into B chains which tilt to accommodate the metal 
atoms (Fig. 3 in Ref. [3ll ). This leads to a more pro- 
nounced relative deviation from Vegard's law^S in oP12- 
FeB 2 (-16.4%) compared to that in hP6-ReB 2 (-6.4%). A 
qualitative analysis of the materials hardness can be done 
by comparing shear moduli, as they have been shown to 
correlate with hardness in some cases^ The shear mod- 
uli in the predicted B-rich phases are at least as high as 
in the known oP8-FeB and tI12-Fe2B compounds (Ta- 
ble I). Although elastic constants alone are not sufficient 
for a quantitative prediction of materials hardness^i it 
is encouraging to find the elastic moduli and the Poisson 
ratio in oP12-FcB2 to be comparable to those in hP6- 



ReB2 (see Tables I and II in Ref. Hi] for comparison of 
experimental and theoretical values). 

In summary, our high-throughput simulations span- 
ning a large library of known and proposed structures 
predict consistently that new compounds should form 
under accessible synthesis conditions in such a common 
binary system as Fe-B. Ground state search under am- 
bient or elevated pressures can be expanded further by 
performing unrestricted structural optimization driven, 
e.g., by evolutionary algorithms^ for larger unit cells and 
other Fe-B compositions. Due to the known exceptional 
complexity of metal boride structures^ experimental in- 
put may be key for determination of the true ground 
states. 
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FIG. 1: (Color online) Calculated formation enthalphy of Fe- 
B compounds at P = 20 GPa. The cyan and red labels cor- 
respond to known and proposed ground states, respectively. 
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FIG. 2: (Color online) Comparison of competing oIlO- and 
oPlO-CrB.4 structures with relaxed lattice vectors of a — 4.749 
A, b = 5.488 A, c = 2.852 A and a = 4.723 A, b = 5.474 
A, c = 2.851 A, respectively: a) simulated x-ray patterns 
for A = 1.5418 A; b) calculated phonon dispersons showing 
dynamic instability of oI10-CrB4 at q — T. 



Library of structure types. We have considered over 
40 commonly seen M-B and M-C structure types listed 
in the Inorganic Crystal Structure Database (ICSD). 1 
The list below is ordered by composition and specifies 
the Pearson symbol and the prototype for each entry. 



1:12: cF52-UBi 2 ; 1:7: oI64-MgB 7 ; 3:20: oS46-Na 3 B 20 ; 
1:6: cP7-CaB 6 ; 3:14: tI160-Li 3 Bi 4 ; 1:4: tP20-UB 4 , 
oI10-CrB 4 , mS10-MnB 4 ; 1:3: tP20-LiB 3 , hP16-Mo . 8 B 3 ; 
2:5: hR21-Mo 2 B 4 . 65 ; 1:2: hP3-AlB 2 , oP6-RuB 2 , hR18- 
MoB 2 , oP12-PbCl 2 , hP6-ReB 2 , hP6-MoS 2 ; 2:3: oS20- 
V 2 B 3 , hP10-Ru 2 B 3 ; 3:4: oI14-Ta 3 B 4 ; 5:6: oS22-V 5 B 6 ; 
1:1: oP8-FeB-b, oS8-PtB . 67) tI16-MoB, 0S8-TII, 0P8- 
TiSi, hP2-WC; 5:4: hP18-Rh 5 B 4 ; 11:8: oP38-Ru n B 8 ; 
3:2: tP10-U 3 Si 2 , aP30-Rh 9 B 5 . 5 ; 5:3: tI32-Cr 5 B 3 ; 2:1: 
tI12-CuAl 2 ; 7:3: hP20-Th 7 Fe 3 , oP40-Mn 7 C 3 ; 5:2: mS28- 
Mn 5 C 2 ; 3:1: oP16-Fe 3 C, oS16-Re 3 B, tI32-Ni 3 P; 23:6: 
cF116-Cr 23 C 6 ; 4:1: tP10-Be 4 B. The structures proposed 
in our previous study 2 - include 1:4: oP10-FeB 4 , 1:3: mP8- 
FeB 3 , 1:2: oP12-FeB 2 . 
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FIG. 3: (Color online) Nearest neighbor histograms in fully 
relaxed boron-rich Re-B and Fe-B structures. In oP10-FeB4 
the neighbors are shown only for one of the two 4g boron sites. 



Comparison to published data. We use available in- 
formation on stable iron and rhenium borides at ambi- 
ent pressure to benchmark the performance of the gen- 
eralized gradient approximation 3 employed in this study. 
For tI12-Fe 2 B and oP8-FeB our calculated magnetic mo- 
ments, 1.85 and 1.16 /x^/Fe, are in good agreement 
with previously measured values of 1.91 and 1.03 //s/Fe, 
respectively^ As shown in Table HI for tI12-Fe 2 B the cal- 
culated heats of formation are close to the previously 
calculated values but are noticeably larger in magnitude 
compared to the measured ones; for oP8-FeB we observe 
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TABLE I: Heat of formation for two stable compounds in 
the Fe-B system ( a) reference states: a-Fe, /3-B, T=298K; 
^ reference states: 7-Fe, /3-B, T=1385K; c ' values originally 
given in eV/atom). 



an excellent agreement between theory and experiment. 
The elastic constants are determined by fitting the total 
energy of distorted unit cells to polynomials of second- 
order in strain. In order to determine all independent 
elastic constants for a given cell type, the same number 
of different strains is used. For each applied strain in 
the range [—0.03 : 0.03] with 0.01 steps the internal co- 
ordinates were fully relaxed. We find a generally good 
agreement between the present and previously reported 
values for the elastic constants and moduli, see Tabic [TTJ 



Cn Cr 



C13 C22 C2 



C: 



33 



C44 C55 Cr 



tI12-Fe 2 B 



theo. 
theo. 
theo. 
exp. 
exp. 
exp. 
exp. 



389.8 169.6 218.5 



413 154 132 



282.2 114.0 



389 148 



157 



B 
249.7") 

331.04 5) 
228 



164 



G 
60.2 a) 
152.77 

143 



E 
184.4 
397.22 
355 
343 
290 
300 



0.388 
0.3 
0.24 



Ref. 
9 
6 

This work 

M 

13 
14 



oP8-FeB 



theo. 
exp. 
exp. 
exp. 
exp. 



371 250 



431 209 505 207 118 194 



287 



138 



358 
284 
350 
600 
125-624 



0.29 



This work 
10 
11 
13 
15 



hP6-ReB 2 theo. 641 159 128 
theo. 631.3 158.4 133.8 



theo. 
theo. 
exp. 
exp. 
exp. 



689 
632 



170 

154 



141 
131 



1037 271 
1015.0 257.0 
1090 281 
1005 248 



350 c) 
347. 7 c) 
369 a) 
339 a) 
360 
382 



283 c) 
273.5 c) 

a) 
a) 



294 
266 



273 



696 
632 



0.19 
0.19 



371±22 302±18 712±43 



16 

XL 

18 

This work 
19 
20 
21 



TABLE II: Comparison of calculated elastic constants and moduli to previously reported theoretical and experimental values 
for selected known compounds. The elastic moduli were calculated using a ' Voigt-Reuss-Hill theorem, 6 - ) equation of state, or 
c ' Voigt theorem. 
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